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ABSTRACT 

We have used the Very Long Basehne Array (VLBA) at 43, 23 and 15 GHz 
to measure the solar gravitational deflection of radio waves among four radio 
sources during an 18-day period in October 2005. Using phase-referenced radio 
interferometry to flt the measured phase delay to the propagation equation of 
the parameterized post-Newtonian (PPN) formalism, we have determined the 
deflection parameter 7 = 0.9998 ± 0.0003 (68% confldence level), in agreement 
with General Relativity. The results come mainly from 43 GHz observations 
where the refraction effects of the solar corona were negligible beyond 3 degrees 
from the sun. The purpose of this experiment is three-fold: to improve on the 
previous results in the gravitational bending experiments near the solar limb; to 
examine and evaluate the accuracy limits of terrestrial VLBI techniques; and to 
determine the prospects and outcomes of future experiments. Our conclusion is 
that a series of improved designed experiments with the VLBA could increase 
the presented accuracy by at least a factor of 4. 
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Subject headings: gravitation — quasars: individual (3C279) — relativity techniques: 
interferometric 



Introduction 



The parame terized post-Newtonian (PPN) formalism was designed by Nordtvedt and 
Will (jWilll I1993I ) to analyze plausible deviations in the general theory of relativity in the 
case of a spherically-symmetric gravitational field of a massive body. The PPN formalism 
introduces several phenomenological parameters that are targets for experimental measure- 
ments, and the propagation of electromagnetic waves is affected strongly by one of them, 7. 
This parameter empirically measures the devia tion frorn the li near approximation of General 
Relativity in which case the value of 7 = 1.0 (lEinsteinlllQlGl ): The New tonia n limit in this 
linear approximation would imply a value 7 = (see equation 7.22 in [Will (119931 )). The 
parameter 7 can be determined precisely and unambiguously by measuring the deflection of 
electro- magnetic radiation by the gravitational field of the sun, with th e GR value of 1.750 ^^ 



at the solar limb. Such experiments have been performed since 1919 (jPyson et al. Ill920l ) 



and the most accurate interferometric observatioi is have shown that 7 is consistent with 



1976 



unity to an accuracy of less than one part in 1000 ( Shapiro 1964 ; Fomalont fc Sramekl 
Robertson et al. I Il99ll : iLebach et al. I Il995l : IShapiro et all |200J) • Recently, iBertotti et al. 



(I2OO3I ) have made use of the Cassini tracking data in 2002 to set the upper hmit of 7 to a 
few parts in 100,000. Somewh at less accuracy has been obtained with another analysis of 



the data (Anderson et a. 



2004), and there are some concerns about the effect of the motion 



of the solar barycenter (IKopeikin et al.l 120071 ). Hence, additional accurate and independent 



measurements are important to determine the value of the fundamental parameter, 7, using 
the radio-interferometric and other techniques. 

With the demo nstrated accuracy of the VLBA to measure relative positions of radio 
sources to 0.01 mas (IFomalont &: Kopeikinll2003l : iBrunthaler et al. II2006I ). the gravitational 
bending could be measured potentially with the radio-interferometric technique to a few 
parts in 100,000, although the coronal refraction when observing within a few degrees of 
the sun produces large path-length changes. In this paper, we present the results from the 
measurement of 7 performed with the Very Long Baseline Array (VLBA) in October 2005, 
and we suggest how more accurate measurements of 7 can be obtained. 



- 3 - 



2. The Observations and Data Reduction 

2.1. The Observational Strategy and Parameters 

Every year in early October, the sun passes in front of the strong radio source 3C279. 
By measuring the change of its angular position, caused by the gravitational bending effect, 
with respect to other nearby sources over a period of several weeks, the gravitational pa- 
rameter 7 can be accurately determined, and we used the opportunity in 2005 to conduct 
a dedicated VLB A experiment. In contrast to several previous interferometric experiments 
that measured the relat ive position of 3C279 with respect to 3C273 using the group delay 



(e.g. lLebach et al. Ill995l ). our observations used the phase delay, a more accurate astrometric 
quantity, to measure the change of relative position of 3C279 with respect to that of several 
fainter sources within a few degrees in the sky. 

For phase referencing observations among close sources, the choice of sources to observe 
must be balanced between two competing factors: the closer th e radio sources a re in the sky. 



the more accurate their relative positions can be determined (iFomalontI l2005l ) : the farther 
apart the radio sources are in the sky, the greater the relative gravitational deflection among 
them. The angular scale pertinent to this balance is related to the closest that successful 
source observations can be made near the sun, about two degrees. The sensitivity of the 
VLBA determines whether there are any sources within a few degrees of 3C279 that can be 
detected in order to measure their accurate posit ions. From an existing catalog of compact 



and relatively bright sources (jPetrov et al. 1 120061 ) . several candidates were found, and after 



brief VLBA observations, three were chosen for this experiment, and their configuration in 
the sky is shown in Fig. 1. The apriori positions and total flux densities of the sources are 
shown in Table 1, and their positions, estimated uncertainties and compactness have been 
determined within the International Celestial Reference Frame system to an accuracy of < 1 
mas. Even the weakest of the sources can be detected by the VLBA within about 20 seconds 
of integration time. 

The sources were observed with the VLBA on eight days (sessions) centered around 
the October 8 solar occultation of 3C279: October 1, 5, 6, 7, 9, 10, 11 and 18, for a 
period of six hours each session when the sources were above 10° elevation at most of the 
VLBA antennas. Because the ionized plasma in the solar corona produces a frequency- 
dependent deflection that increases inversely with wavelength squared, observations at the 
highest routinely available VLBA frequency of 43 GHz were chosen to minimize the coronal 
deflection. We also observed the sources at 23 GHz and 15 GHz in order to estimate the 
coronal deflection that was present. 

Since the VLBA cannot observe simultaneously among 15, 23 or 43 GHz, it is not pos- 
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sible to remove the coronal refraction instantaneousljQ. Thus, we had to choose a switching 
cycle among the sources and frequencies which minimized the effects of phase fluctuations 
from the troposphere and corona. Because the short-term troposp heric phase fluctuations a t 
43 GHz can limit the coherence time-scale to less than a minute (IBeasley fc Conwaylll995l ). 
it is crucial to switch observations (scans) among the sources within this period in order to 
keep the phase coherent between source scans. Since a change in the observing frequency 
at the VLBA takes at least 15 sec, any frequency switching would have considerably length- 
ened the time between scans, beyond the nominal coherence time. Thus, we decided switch 
among the sources at one frequency with 40-sec scans for about 20 minutes. Then, change 
frequencies and switch among the sources for another 20 minutes. This observing scheme is 
illustrated in Fig. 2. Essentially, each 20-minute period of source switching can determine 
their accurate relative positions at one frequency. The slower cycling among the frequen- 
cies could determine the different positions among the frequencies caused by the coronal 
refraction, albeit on a relatively slow time scale of about one hour. 



The relevant observational parameters were as follows (jWalkeii Il995l ). At each of the 
three observing bands, we observed at four frequencies (IF's), each with a bandwidth of 16 
MHz: 14.93, 14.98, 15.20, 15.37 GHz; 23.02, 23.07, 23.29, 23.46 GHz; 42.73, 42.78, 43.00, 
43.17 GHz. Each IF was separated into 16 frequency channels, each of width 1 MHz, and 
then sampled with two bits. The separation of the frequencies at each band permitted the 
determination of the group delay (phase slope versus frequency) to improve the tropospheric 
delay estimates using many calibrator observations at the beginning and end of each day. 
We also included several observations of the strong sources 3C273 and J1310+334 over the 
day to monitor the instrumental delay changes. The observing schedule for the eight sessions 
was identical apart from slightly different start and stop times. 



2.2. Initial Data Calibration, Editing and Averaging 

The data from all antenna-pairs were processed using the VLBA correlator in Socorro, 
NM, and the output visibility data were averaged into 2-sec samples. The apriori model in the 



correlator uses the NASA CALC software package (http://www/gemini.gsfc. nasa.gov/solve) 



that includes the gravitational bending, assuming 7 = 1. Hence, the fitting parameter used 



^The 2.3/8.4 GHz simultaneous frequency system, whic h is available on the VLBA and used for most 
astrometric observations to remove the ionosphere refraction ( Ma et al. 1998h . is subject to severe scattering 



by the coronal refraction within about four degrees of the sun. Previous observations to measure the solar 
deflection at these low freque ncies alone did not produce 7 accuracies better than the 10^^ level. Hence, 
higher frequencies were used ( Lebach et al. IQQSh . 
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in the data analysis is the departure of 7 from unity. 



The standard VLBA cahbrations and editing were made on the correlated data (lUlvestad. Greisen &: M 



2OOII ) independently for 15 GHz, 23 GHz and 43 GHz: 1) The apriori amplitude calibrations 
of each antenna, from the monitored system temperature and aperture efficiency, as well as 
the correlator bias, were applied; 2) Data were flagged during receiver/antenna malfunctions, 
antenna mis-pointing, and when the source elevation for any antenna was below 10°. 3) Im- 
proved parameters for the earth-orientation and rotation, polar motion, antenna motions and 
nutation were available from on-line services and applied to the correlated data; 4) The cal- 
ibrator observations at the beginning and end of each session were analyzed to determined 



corrections to the zenith path delay used in the correlator model iMioduszewski Sz Kogan 



(|2005|) 



The data from the correlator for each observing frequency contained 64 streams: four 
IF frequencies, each split into 16 frequency channels. The relative gains and phases among 
the streams change slowly and were determined using several observations in each session 
of strong sources, 3C273 and J1310+332 (and 3C279 when not too close to the sun). The 
relative gain and phase determination among the frequency streams varied about 5% and 
15°, respectively, over a session, so the changes were easily corrected. The 64 data streams 
were then coherently averaged to form one data stream with high signal to noise ratio. These 
data at 15 GHz, 23 GHz and 43 GHz were then used to determine the shorter time-scale 
phase changes induced by the troposphere and solar corona, described in the next section. 



2.3. Tropospheric Editing 

When observing at high radio frequencies, small "clouds" of water vapor pass through 
each antenna beam and produce variable delays of a few millimeters (about 100° at 43 GHz) 
over time scales of a minute or less. These induced phase variations limit the astrometric and 
imaging accuracies. An example of relatively stable phase behavior is shown in Fig. 3a. Since 
3C279 is a strong source, its phase can be determined for each 2-sec sampled data point over 
a 40-sec scan (closely-spaced points), and over this period the typical phase variation over a 
scan is about a 30°. In order to have co/ierent observations, the phase between the consecutive 
3C279 scans must be unambiguously connected with no lobe ambiguities, otherwise the 
residual phase of the other sources(their measured phase minus the interpolated phase from 
3C279) will be grossly in error. For the BR and MK baselines, the phase connection between 
scans of 3C279 is obvious, although residual phase offsets of 50° do occur. However, for the 
SC baseline, the scan to scan phase changes of 3C279 are larger and in some cases (between 
18:41 to 18:42) it is not clear whether the phase increased or decreased between scans, so 
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this period with SC is of dubious astrometric quahty. 

Using plots similar to that of Fig 3a for all observations and frequencies, we edited 
periods when the phase coherence was poor. These periods often were associated with 
inclement weather conditions at an antenna site. The percentage of data removed was 4% 
at 15 GHz, 7% at 23 GHz and 17% at 43 GHz, mostly for the Saint Croix VI and Hancock, 
NH antennas that are at relatively humid locations. 



2.4. Coronal Editing 

In Fig. 3b we show a 2-min segment of data at 23 GHz for 3C279 on October 9 when the 
source was about 1.7° from the sun. Although the phases are coherent over 2 seconds, there 
are rapid changes of phase over 10 seconds or longer, hence the phase cannot be connected 
between scans. Notice that the phase changes for the shortest baseline (bottom plot) of 
240 km are as variable as those for the longer baselines. This indicates that the coronal 
refraction significantly varies over second time-scales with a linear scale-size smaller than a 
few hundred kilometers in the solar coron£@. 

We found that when a source was closer than ~ 3° from the sun, phase-stable observa- 
tions often were not possible because of the severe and short term coronal refraction at all 
frequencies. Thus, all observations on October 6 and 7, when three of the four sources were 
within this angle from the sun, were astrometrically worthless. Other data not used because 
of coronal turbulence were for J1248 and 3C279 on October 5 and 3C279 on October 9. 



3. Astrometric Analysis 

3.1. The Removal of the Source Structure 

Although all of the sources are dominated by a compact core component, extended 
emission is also present and often associated with a jet component that emanates from the 
core. This structure causes two astrometric problems: the peak brightness of the core (best 
definition of the location of the source) is resolution dependent, and the structure may vary 
with time. We, thus, determined the source structures at each frequency on October 1 



^the analysis of the statistical properties of the phase data when the sources were close to the sun can 
provide information on the velocity, density and size of the coronal turbulent cells. This will be reported 
elsewhere. 
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and October 18 using the self-calibration algorithm (ICornwell &: Fomalont I Il989l ). Since 
this algorithm gives no positional information about the source, we arbitrarily placed the 
maximum brightness for each source at each frequency at the assumed apriori position of 
the source given in Table 1. This convention has no effect on the astrometric results since 
the analysis determines the residual true sky position of the peak of the source, relative to 
the assumed model and apriori position. 

Because significant structure and intensity changes in most extra-galactic source radio 
emission occur o ver months of time , little change was expected over the 17-day period of 



this experiment (iHughes et allll992l ). Changes in structure can also be indirectly inferred 
also by the variability of the peak flux density of the sources that are given in Table 2. The 
only significant variability occurred for 3C279 and other monitoring showed that the source 
reached a maximum of 20 Jy at 23 GHz in early 2005 and decreased through the year. Hence, 
the decrease in flux density is real and the structure did change slightly over the experiment. 
Comparison of the accurate structures for 3C279 on October 1 and October 18 showed that 
minor variations occurred near the position 0.2 mas west and 0.1 mas east of the radio core 
at 43 GHz, but that the apparent position change of the peak brightness at 43 GHz and 23 
GHz was less than 0.02 mas. 

The source structures will be discussed and displayed in more detail elsewhere (Fomalont 
et al, in preparation), but a summary follows: The sources J1246 and J1248 are nearly point 
sources at 43 GHz, with a faint component at 15 GHz about 0.5 mas to the west and 1.4 
mas to the north-east, respectively. The source J1304 is a nearly equal double at 43 GHz 
with a separation of 0.8 mas, and 3C279 has a jet component which extends about 1 mas 
south-east of the core. 

Since the source structures were essentially unchanged over the 18-day period of the 
experiment, we removed their effects by dividing the visibility data with the visibility data 
model associated with the source structure (averaged from the Oct 1 and Oct 18 images) at 
each frequency. This produced a structure-free (effectively now a point source) data set for 
each source and frequency. Any resolution differences among the observing sessions will not 
produce a change in the location of the peak brightness of this revised data set. 



■^http:/ /www. via. nrao.edu/astro/calib/polar/2005 
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3.2. Phase Referencing to Obtain the Residual Phase 

After the initial calibrations, editing, and removal of the source structure, we then used 
the standard phase referen cing technique on the v isibility data to determine the relative 



position among the sources (IBeasley fc Conwaylll995l ). This technique uses one of the sources 



as the main reference to determine the phase error associated with each antenna for each 
reference source scan. This phase error is then interpolated between each reference scan, and 
applied to the other sources to obtain their residual phases. This technique is illustrated in 
Fig. 3a. Since the switching time and proximity in the sky among the sources were small, 
the phase errors associated with the reference source are virtually the same as that for all 
sources. Hence, any residual phases found for the non-calibrator sources are predominantly 
associated with a position offset from its apriori value. We used 3C279 as the main phase 
reference on October 1, 11 and 18 when it was not too close to the sun, and J1304 as the 
main phase reference on the other days. The phase referencing was done independently for 
each frequency. 



3.3. Source Position Determination from the Residual Phase 

Two analysis methods were used to determine the radio source positions from the mea- 
sured residual phases for each source and frequency. The first method used the basic in- 
terferometric Fourier imaging and deconvolution techniques to produce an image of each 
source. We then determined the location of the peak brightness of the source (all nearly 
point sources since the structure component has been removed) , and an error estimate based 
on the quality of the image and the signal-to-noise ratio at the peak intensity. For the second 
method, we fit the measured residual phases directly to a point source model and determined 
by a least-square analysis the position offset from the apriori position and the error estimate 
based on the phase deviations from the best fit. 

Although the relative source positions at each frequency could be derived from every 
20-min segment of data, the expected change of position during each session caused by the 
coronal bending and by any small departures of 7 from unity was estimated to be less than 
0.20 mas at 15 GHz, and smaller at the higher frequencies. Hence, the data at each of the 
frequencies for an entire session could be analyzed with no significant position smearing. 
However, for 3C279 on October 10 when the coronal bending significantly changed, the 
analysis was split into two half-session periods. 

Both methods used the same set of data, and they provided somewhat different positions 
and errors that are associated with the different effective weighting of the data points. For 
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the source position, we used the average position obtained from the two methods. The 
methods provided three position errors, one from each analysis method and one also from 
the difference between the two position determinations. We chose the largest of the three as 
the error estimate. Finally, we adopted a floor uncertainty that was based on the theoretical 
signal-to-noise ratio of the data: 0.02 mas E/W and 0.03 mas N/S for 3C279, and 0.05 ms 
E/W and 0.07 mas N/S for the other three sources. 

The resulting source positions from the above analyses, made independently for each 
frequency, are shown in Table 3. In order to have a one fiducial reference point over the 
entire experiment, we have listed in column 1 the relative positions of 3C279, J1254 and 
J1256 with respect to J1304 for all days. This requires only a trivial position translation to 
J1304 of the positions on days when 3C279 was used as the phase referenc^. The distance 
of the source from the sun in degrees in mid-session is shown in column 2. When one or 
more of the sources were within 2.5° of the sun, the coronal turbulence was often too severe 
to obtain phase coherence and no positions could be determined. For this reason, there are 
no October 6 and 7 entries in the table. On October 10 when 3C279 was relatively close to 
the sun, the analysis was made on the first half and second half of the session because of the 
relative large change in coronal refraction on this day. 

For observations on October 1 and 18, when the sources were sufficiently far from the 
sun, the maximum expected relative gravitational bending was 0.02 mas and coronal bending 
at 15 GHz was 0.03 mas, less than the position errors for each day. Hence, the positions 
on these two dates were chosen to be the undefiected relative positions at each frequency 
among the sources. The position differences between the two days are generally consistent 
with the estimated errors, although there are several outliers are about 3-cr, especially for 
the N/S positions for J1246 and J1248. 

3.4. Position Changes with Session and Frequency 

Table 4 lists the position changes from the non-defiected positions for the three sources 
with respect to J1304, as a function of session and frequency. The assumed non-deflected 
position is given on the first line for each source and frequency, and it is the weighted 
average of the October 01 and October 18 positions, given in Table 3. The offsets of about 



^Because the experimental errors are dominated by antenna-based residual tropospheric and ionospheric 
errors, rather than by receiver noise and other baseline-stochastic processes, there is no additional information 
by using the results from all six source-pairs, compared with the non-degenerate set of the three source-pairs 
listed. 
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1 mas for 3C279, J1246 and J1248 reflect the error in the apriori position of J1304. The 
offset positions among 3C279, J1246 and J1248 are in good agreement and suggests that their 
apriori positions are consistent to about 0.3 mas. The shght differences in the offset positions 
among the three frequencies are real and show the affect of sturcture changes for each source 
among the frequencies. A more complete discussion of the registration of the source positions 
among the three frequencies is given elselwhere (Fomalont et al, in preparation). 

Below this line, the table then lists the position changes from the undeflected positions 
for the other observing dates. Since the gravitational bending with 7 = 1 between the 
sources was included in the correlator model, the position changes will include only the 
gravitational bending associated with a (7 — 1) term, as well as that caused by noise and 
remaining tropospheric and coronal refraction components that differ among the sources. 

The first line for an observing date, labeled Each Freq, lists the residual positions of 
the sources for each of the three frequencies. All entries are obtained from the relevant 
entries in Table 3. The second line for an observing date, labeled Corona Free, gives the 
residual coronal-free position. Since the c oronal plasma deflectio n varies with the observation 



wavelength squared at radio frequencies (IThompson et alll200ll ). the estimate of the coronal 



free position, po, can be obtained from the relationship. 



Pii^i) = Po + ll^i 



where p(t'j) are the measured positions (x,y) at the three frequencies z/j, and I is the average 
coronal refrac tion (magnitude and direction in the sky) at 1 GHz for the observation period 
( Rogersl 1970 ). We determined the estimated values and uncertainty of the ionosphere cor- 
rection and the corona-free position of the sources using a least-square fit. The discussion of 
the effectiveness in obtaining corona-removed solutions is given in the next section. 



4. The Determination of the Parameter 7 

The approximate accuracy for the determination of 7 from Table 4 can be estimated in 
a straight-forward manner. Since the typical position uncertainty of one measurement (one 
source-pair per day) is ~ 0.07 mas for the differential gravitational bending of ~ 150 ma^, 
the measurement accuracy of 1 part in 2000 corresponds to an accuracy of 7 of 1 part in 



^For example, On October 11 the average relative deflection between 3C279 and J1304 was 100 mas and 
66 mas in the e/w and n/s directions, respectively 
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1000. Since about 20 independent measurements of the bending were made (10 source pairs 
with two position coordinates), a sensitivity of 7 to a few parts in 10~^ is expected. 

The more formal determination of the optimum value of (7 — 1) from the experimental 
results in Table 4 is straight-forward. The analysis minimize the (normalized) chi-squared 
expression 

2 f Pd{^)-0.5h-l)Dd(^) Y 
arf(z) J 

d,i 

where Pd{i) and ad{i) are the measured position offset and error estimate, respectively, from 
Table 4. The {i) loops over the sources (3C279, J1246, J1248), and the (d) loops over each 
day or half-day observation (Oct05, 09, 10a, 10b, 11). The term -Drf(i) is the differential 
General Relativity gravitational bending prediction, averaged over the session. The sum 
is made over the E/W and N/S values and estimated errors separately, so the number of 
degrees of freedom is k = 19, if all ten observing points are used. 

The results for several solutions for 7 using the different data sets from Table 4 are listed 
in Table 5. The normalized indicates the ratio of the rms of the best fit divided by that 
expected from the estimate error of each entry. The 43 GHz corona-free fit has the lowest 
for two reasons: the lessening of some coronal effects, and the increase of the position 
errors. The two 43-GHz only solutions (with no removal of the ionosphere contribution) 
show the effect of the Oct05 session that was made relatively close to the sun. Finally, the 
23-GHz only solution has a relatively large normalized x^, even excluding the Oct05 session, 
and suggests that coronal refraction, which is four times larger than that at 43 GHz, is 
dominating the sensitivity of the experiment at 23 GHz. Nevertheless, the variation of 7 and 
its estimated error among the different solutions are in good agreement and suggest that the 
determination of 7 change little between several different analyses of the data. 



5. Results and Discussion 



The solutions in Table 5, using different selections of data, are all consistent with GR 
value of the parameter 7=1. For the results in this paper, we have taken an average of 
the four solutions to obtain 7 = 0.9998 ± 0.0003. This result is a factor of three times more 
accurate than the previous dedicated radio interfer ometric observations made specifically to 
measure the gravitatiori al bending: 0.9996± 0.0017 Jlebach et al. Ill99,'^ . and 1.0002±0.0010 
( jRobertson et al. Ill99ll ). Usi ng the radio geodetic data base of thousands of observing ses- 
sions between 1979 and 1999, Ishapiro et al j^OoJ obtained 7 = 0.9998 ± 0.0004, compara- 
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ble in accuracy to the present results. The result from the measurem ents as the spacecraf t 
Cassini passed by the sun in September 2002 is 7 = 1.00002 ±0.00002 JSertotti et allboosh . 
and was discussed in §1. 

Although expected departures from 7 = 1 are likely to be a factor of 10 to 100 smaller 
than the accuracy of this experiment, some possibilities are: (1) The long-range scalar 



tensor interaction in scalar-tensor theories of gravity (j Willi 1 19931 ). and lDamour fc Nordtvedt 



(I1993I ) predict a lower bound for the present value of 7 at the level of 10 ^ ~ 10 ^; (2) 



The long-range vector-tensor interaction in vector-tensor theories of gravity (IKostelecky 
2004 ) may suggest a "spontaneous vio lation" of the Lorentz invariance and could modify the 
value of 7 (IBailey fc KosteleckyI 120061 ) : (3) The more complicated nature of the gravitational 
coupl ing between the curvatu re and stress-energy tensor of matter may lead to changes 
of ^ (I Jaekel fc RevnaudI 120051) : (4) The plausible existence of the effective graviton's mass 
teabak fc Grishchukll2003h i ;hat would avoid the van Dam-Veltmann-Zakharov discontinuity 
(Ivan Dam fc Veltmarull97d : IZakharovlll970l ) would also effect the value of 7. 



There are several changes in the design of a similar VLBA experiment that should 
improve the accuracy of 7 by about a factor of two. First, by choosing a set of sources 
that can be observed when the sun is further nortlfl, each day's integration time can be 
increased from 6 to 10 hours, and the sources will be on average at higher elevations than 
those from the present experiment. We estimate that a more northern experiment will lower 
the position rms by about 20%. 

Secondly, we lost some sensitivity in the determination of 7 by not observing on days 
when the gravitational deflection was relatively large (egs October 3, 4, 12, 13). We found 
that the most accurate results are obtained between solar distances of three to five degrees, 
corresponding to about three days to seven days from the nearest approach of the sources 
with the sun. Thus, the addition of the above four observing days in the 2005 experiment 
would have decreased the position rms by 25%. 

Finally, the 2005 experiment devoted too much observation time at 15 GHz and possibly 
23 GHz in an attempt to lessen the effect of coronal bending. When the coronal bending 
was significant at 43 GHz (> 0.1 mas), the phase stability degraded considerably at all 
frequencies and no astrometric information could be obtained. Thus, future observations 
should concentrate on 43 GHz, somewhat further from the sun. Additional observations at 
23 GHz are useful to monitor the coronal refraction since it can be variable with time, but 
observations at 15 GHz are not helpful in determining the coronal refraction. The potential 



^Several possible experiments are: March 30 to April 24— J0121+1149, J0104+1134, J0139+0842, 
J0129+1146; April 10 to May 01— J0204+15114, J0158+1307, J0209+1352, J0211+1051 
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additional observing time by a factor of two at the primary frequency of 43 GHz in a new 
experiment, compared with the 2005 experiment, would decrease the position rms by about 
30%. 

Hence, the above three changes of strategy in a VLBA-designed experiment will clearly 
improve the accuracy of the determination of 7 from a single 10-day experiment by at least 
about a factor of two compared with that of the October 2005 experiment. The availability 
of many such experiments over the year, all repeatable from year to year, can provide in- 
dependent results estimates of 7 since the dominant error is produced by the quasi-random 
errors of the short-term troposphere and, to a lesser extent, coronal refraction. Also, any 
small systematic errors associated with the different sky configuration for the sources among 
the other possible experiments should be somewhat independent. Thus, we expected that 
the accuracy of 7 should improve by roughly the square-root of the number of experiments. 
Hence, with sufficient VLBA resources, the uncertainty in the parameter 7 can be decreased 
by at least a factor of four compared with the results given in this paper. 
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Fig. 1. — The Source Configuration for the Deflection Experiment: The solar tra- 
jectory between October 1 and 18 is shown by the diagonal line, with the eight observing 
days superimposed. The location of the four radio sources are indicated. 
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Fig. 2. — The Observing Strategy: The top plot shows the switching among the four 
sources within every 20-minute period. The bottom plot shows the switching of these 20-min 
blocks among the three frequencies through each observing day. 
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(a) Troposplicric Phases on October 1 



(b) Coronal Phases on October 9 
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Fig. 3. — (a) Tropospheric Phase Behavior: The phase behavior at 23 GHz on October 
1, when the sources were > 5° from the sun, is shown for three VLB A baselines from Los 
Alamos, NM (LA) to: Brewster, WA (BR), 1800 km; Mauna Kea, HI (MK), 5000 km; 
and Saint Croix (SC), VI, 4500 km. The period covers 13 minutes near the middle of the 
day, with the UT time given at the bottom. The densely-packed points are the phases 
for 3C279 at two-second intervals, and the more isolated larger points show the phases for 
the other calibrators, averaged over 20 sec (two per scan). All sources he on a relatively 
continuous temporal phase curve because their accurate relative positions were determined 
before producing the plot, (b) The Coronal Phase Behavior: The phase behavior at 
23 GHz on October 9 for 3C279, when the source was 1.7° from the sun, is shown for three 
VLBA baselines over a 2-min period. The baseline between LA and Pic Town, NM (PT) 
is only 240 km. Each plotted point shows the measured phase for 3C279 at two-second 
intervals. The other signal from the other three sources observed during the gaps were too 
decorrelated to be detected. 
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Table 1. Apriori Radio Source Position and Total Flux Density 



Source 


RA DEC 


S43GHZ 


S23GHZ 


SwGHz 




Equinox 2000.0 




Jy 




3C279 


12 56 11.166557± 0.000013 -05 47 21.52467± 0.00031 


11.22 


12.75 


15.27 


J1246-0730 


12 46 04.232100± 0.000014 -07 30 46.57456± 0.00031 


0.24 


0.34 


0.50 


J1248-0632 


12 48 22.975657± 0.000016 -06 32 09.81737± 0.00041 


0.16 


0.19 


0.26 


J1304-0346 


13 04 43.642255± 0.000022 -03 46 02.55122± 0.00065 


0.38 


0.48 


0.65 
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Table 2. Peak Flux Density Variations of Sources 





Date/Freq 




T1 9/lfi 


J1248 




J1304 






Oct 01:43 GHz 


lO.li 0.3 


0.50± 0.02 


0.40± 0. 


,02 


0. 


,44± 


0. 


,03 




Oct 18 


:43 GHz 


9.0± 0.4 


0.54± 0.02 


0.48± 0. 


,02 


0. 


,51± 


0. 


,03 




Oct 01 


:23 GHz 


12.9± 0.2 


0.61± 0.02 


0.41± 0. 


,03 


0, 


,52± 


0. 


02 




Oct 18 


:23 GHz 


11.2± 0.4 


0.63± 0.02 


0.50± 0. 


,04 


0. 


,57± 


0. 


03 




Oct 01 


:15 GHz 


13. 7± 0.7 


0.62± 0.02 


0.28± 0. 


,02 


0. 


,51± 


0. 


03 




Oct 18 


:15 GHz 


12.9± 0.3 


0.59± 0.02 


0.27± 0. 


,02 


0. 


,49± 


0. 


,04 






Table 3. 


Source Positions with Respect to J1304 






Source 




43 GHz 


23 GHz 






15 GHz 


Date 


(deg) 


E/W 


N/S 


E/W 


N/S 




E/W 




N/S 


3C279 






















OctOl 


6.6 


0.96 ± 0.03 


-0.11 ±0.04 


0.77 ± 0.03 


0.16 ±0.04 




0.78 ± 0.04 




-0.05 ± 0.05 


Oct05 


2.8 




















Oct09 


1.3 




















OctlOa 


2.2 


0.95 ±0.10 


-0.09 ±0.14 


0.65 ±0.14 


-1.27 ±0.20 




0.36 ±0.23 




-2.30 ± 0.30 


OctlOb 


2.4 


0.93 ±0.13 


-0.17 ±0.18 


0.56 ± 0.16 


-1.28 ± 0.24 




0.20 ±0.28 




-2.80 ±0.40 


Oct 11 


3.2 


1.03 ±0.03 


-0.11 ±0.04 


1.02 ± 0.04 


-0.05 ±0.05 




1.11 ±0.07 




-0.29 ±0.09 


Octl8 


10.2 


0.93 ± 0.02 


-0.08 ± 0.03 


0.83 ± 0.02 


0.16 ±0.04 




0.84 ± 0.04 




-0.02 ± 0.05 


J1246 






















OctOl 


5.5 


0.75 ± 0.05 


0.16 ±0.07 


0.64 ± 0.05 


0.54 ± 0.07 




0.81 ±0.05 




0.33 ± 0.07 


Oct05 


2.6 


0.46 ± 0.06 


0.15 ±0.10 


-0.23 ± 0.08 


0.49 ±0.11 




0.22 ±0.10 




2.15 ±0.14 


Oct09 


3.7 


0.76 ± 0.05 


0.09 ± 0.07 


0.82 ±0.07 


0.28 ±0.11 




0.78 ± 0.05 




0.16 ±0.07 


OctlO 


4.5 


0.76 ±0.05 


0.09 ± 0.07 


0.82 ± 0.07 


0.28 ± 0.11 




0.78 ± 0.05 




0.16 ±0.07 


Octll 


5.4 


0.75 ± 0.08 


-0.09 ±0.12 


0.66 ±0.10 


0.29 ±0.14 




0.75 ±0.15 




0.18 ±0.10 


OctlS 


12.1 


0.71 ± 0.05 


0.18 ±0.07 


0.68 ± 0.05 


0.38 ± 0.07 




0.86 ± 0.05 




0.21 ± 0.07 


J1248 






















OctOl 


5.3 


0.70 ± 0.05 


-0.40 ±0.07 


0.66 ±0.05 


-0.02 ±0.07 




0.85 ± 0.05 




-0.15 ±0.07 


Oct05 


1.8 




















Oct09 


3.0 


0.68 ±0.09 


-0.32 ±0.11 


0.80 ± 0.18 


0.00 ± 0.23 




1.02 ±0.10 




0.10 ±0.14 


OctlO 


3.9 


0.71 ±0.13 


-0.13 ±0.15 


0.65 ± 0.13 


-0.19 ± 0.19 




0.59 ±0.21 




1.19 ±0.33 


Octll 


4.9 


0.70 ± 0.08 


-0.49 ±0.12 


0.74 ±0.10 


-0.40 ±0.15 




0.64 ±0.15 




-0.84 ±0.20 


OctlS 


11.8 


0.65 ± 0.05 


-0.41 ±0.10 


0.67 ± 0.05 


-0.25 ±0.07 




0.85 ± 0.05 




-0.21 ± 0.07 
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Table 4. Deflection Changes Between the Sources and J 1304 



DATE 


Type 


43 GHz 


23 GHz 


15 GH 


z 






E/W 


N/S 


E/W 


N/S 


E/W 


N/S 






(mas) 


(mas) 


(mas) 




3C279 
















nondcflcctcd position 


0.94 


-0.09 


0.80 


0.16 


0.81 


-0.04 


OctlOa 


Each Freq 


0.00 ±0.10 


0.00 ±0.14 


-0.15 ±0.14 


-1.43 ±0.20 


— U.45 ± U.zo - 


o nil _l_ n Qn 
-z.zb ± U.3U 


Oct 10a 


Corona-Pree 


0.06 ±0.11 


0.39 ± 0.20 










OctlOb 


Each Preq 


-0.01 ±0.13 


-0.08 ±0.18 


-0.24 ±0.16 


-1.44 ±0.24 


-0.61 ±0.28 - 


-2.76 ± 0.40 


OctlOb 


Corona-Pree 


0.06 ±0.14 


0.33 ± 0.25 










Octll 


Each Preq 


0.08 ± 0.03 


-0.02 ± 0.04 


0.22 ± 0.04 


-0.21 ± 0.05 


0.30 ± 0.07 - 


-0.25 ± 0.09 


Oct 11 


Corona-Pree 


0.04 ± 0.05 


0.03 ± 0.04 










J1246 
















nondcflcctcd position 


0.73 


0.17 


0.66 


0.46 


0.84 


0.27 


Oct05 


Each Preq 


-0.27 ±0.06 


-0.02 ±0.10 


-0.89 ±0.08 


0.03 ±0.11 


-1.06 ±0.10 


1.88 ±0.14 


Oct05 


Corona-Pree 


-0.10 ±0.15 


-0.11 ±0.15 










Oct09 


Each Preq 


0.03 ±0.05 


-0.08 ±0.07 


0.16 ±0.07 


-0.18 ±0.11 


-0.06 ± 0.05 - 


-0.11 ±0.07 


Oct09 


Corona-Pree 


0.02 ±0.06 


-0.05 ±0.07 










OctlO 


Eacii Frcq 


0.00 ±0.10 


0.19 ±0.13 


0.06 ± 0.13 


-0.41 ±0.18 


0.11 ±0.18 - 


-0.34 ±0.31 


Oct 10 


Corona-Pree 


-0.02 ±0.11 


0.28 ±0.14 










Octll 


Each Preq 


0.02 ± 0.08 


-0.26 ±0.12 


0.00 ±0.10 


-0.17 ±0.14 


-0.09 ±0.15 - 


-0.09 ±0.10 


Octll 


Corona-Pree 


0.03 ± 0.09 


-0.23 ±0.12 










J 1248 
















nondcflcctcd position 


0.67 


-0.41 


0.67 


-0.14 


0.85 


-0.18 


Oct09 


Each Frcq 


0.01 ±0.09 


0.09 ±0.11 


0.13 ± 0.18 


0.14 ±0.23 


0.17 ±0.10 


0.28 ±0.14 


Oct09 


Corona-Frcc 


-0.02 ±0.11 


0.05 ±0.11 










OctlO 


Each Preq 


0.04 ±0.13 


0.28 ±0.15 


-0.02 ±0.13 


-0.05 ±0.19 


-1.44 ±0.21 


1.37 ±0.33 


OctlO 


Corona-Pree 


0.09 ±0.16 


0.20 ± 0.18 










Octll 


Each Preq 


0.03 ± 0.08 


-0.09 ±0.12 


0.07 ±0.10 


-0.26 ±0.15 


-0.21 ±0.15 - 


-0.66 ±0.20 


Octll 


Corona-Pree 


0.03 ± 0.09 


-0.01 ±0.13 
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Table 5. Solutions for 7 



Solution Type 


7-1 


(7, 


xl 




10-4 


10-4 




43 GHz data (corona-free) 


-2.4 


3.2 


0.9 


43 GHz data only 


-1.0 


2.6 


2.2 


43 GHz data only - Oct5a 


-3.2 


2.8 


1.1 


23 GHz data only - OctSa 


-2.0 


2.4 


4.7 



